Sudden collapse of a colloidal gel 
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Metastable gels formed by weakly attractive colloidal particles display a distinctive two-stage 
time-dependent settling behavior under their own weight. Initially a space-spanning network is 
formed that for a characteristic time, which we define as the lag time rj, resists compaction. This 
solid-like behavior persists only for a limited time. Gels whose age t w is greater than t& yield 
and suddenly collapse. We use a combination of confocal microscopy, rheology and time-lapse video 
imaging to investigate both the process of sudden collapse and its microscopic origin in an refractive- 
index matched emulsion-polymer system. We show that the height h of the gel in the early stages of 
collapse is well described by the surprisingly simple expression, h(r) = ho — At 2 , with ho the initial 
height and r = t w — the time counted from the instant where the gel first yields. We propose that 
this unexpected result arises because the colloidal network progressively builds up internal stress 
as a consequence of localized rearrangement events which leads ultimately to collapse as thermal 
equilibrium is re-established. 

PACS numbers: 82.70.Dd, 83.80.Kn, 64.75.Xc, 83.50.-v 



I. INTRODUCTION 

Soft materials such as colloidal suspensions and emul- 
sions form a remarkably rich variety of non-ergodic states 
[T]-[3] - examples of which are familiar to us in our daily 
life in products as diverse as foodstuffs, surface coatings, 
fabric conditioners, and pesticides. Out-of-equilibrium 
phases occur when suspensions are quenched deep into 
a region of thermodynamic phase separation. Unable 
to phase separate, amorphous solids form which are me- 
chanically rigid but without the long-range translational 
order characteristic of crystalline solids. Slow relaxation 
dynamics prevents the system from reaching their un- 
derlying global equilibrium configurations so these amor- 
phous solids evolve slowly in a complex energy landscape 
with a high number of local minima and as a result dis- 
play glassy dynamics with a rich phenomenology of ef- 
fects such as aging, non-linear responses, and spatial and 
temporal dynamic heterogeneities. 

One of the most dramatic macroscopic manifestations 
of aging is the phenomenon of sudden network collapse in 
gels. Gels consist of a network of particles linked together 
by long-lived attractive bonds. Sedimentation or cream- 
ing of the particles within a gel imposes a buoyant stress 
on the network which since gels are intrinsically rather 
delicate has dramatic consequences for microscopic struc- 
ture and dynamics @HB]. Weak gels, where the strength 
of the attractive potential at contact U c is only a few 
/cbT, show, for instance, a very unusual mechanical re- 
sponse. Initially, the gel behaves as a solid but after a 
finite lag time ra, the gel yields and catastrophically col- 
lapses. Sudden or 'delayed' network collapse is observed 
in a wide variety of materials 7 16 and seems to be ubiq- 
uitous at small U c /k^T. However, while sudden collapse 



* Corresponding author: P.Bartlett@bristol.ac.uk 



has been attributed to channel formation within the gel 
P2J [13] , the microscopic processes operating have never 
been fully established. A better microscopic understand- 
ing of the origin of sudden gel collapse is important not 
only because the distinctive settling behavior is intrigu- 
ing from a scientific viewpoint but also because a quanti- 
tative prediction of gel stability is a critically important 
issue in the formulation and manufacture of many com- 
mercial products. 

The aim of this paper is to report a detailed experi- 
mental study of the stability of gels under gravitational 
stress. We use a colloidal suspension of nearly monodis- 
perse emulsion drops of radius a suspended in an index- 
matched mixture of solvents, which has been well charac- 
terized elsewhere 17]. Gelation is induced by long-range 
attractive depletion forces. Using time-lapse video imag- 
ing we measure the dependence of the height h of a gel 
upon its age t w , counted from the moment when the gel 
was formed. No macroscopic sedimentation is observed 
initially but after a period of latency the gel undergoes 
a rapid collapse as the system separates into colloid-rich 
and colloid-poor phases. We investigate the collapse dy- 
namics as a function of the strength —U c /k^T of the at- 
tractive interactions and the initial height ho of the gel. 
Remarkably we find that when collapse starts the change 
in the height of the gel Ah = ho — h follows a simple 
universal dependence on t w which is independent of the 
initial height ho of the gel. The observation of height- 
independent collapse is surprising and contrasts with the 
marked height dependence seen in short-range gels stud- 
ied to date [T3l HSl [19] . Using a combination of rheol- 
ogy, confocal microscopy and time-lapse video imaging 
we speculate that the collapse of the gel network occurs 
as a result of irreversible aging, via a spatially heteroge- 
neous process of localized 'micro-collapses', which leads 
to a build up of internal stress within the gel and its 
ultimate failure. 

The paper is organized as follows: Section [TT] discusses 



2 



the preparation of the emulsion gels studied and the ex- 
perimental techniques used. Section [Hi] details experi- 



mental results from both macroscopic and microscopic 
measurements on the settling behavior of suspensions of 
attractive particles. The interpretation of the results in 
terms of internal stress relaxations is discussed in Sec.llVl 
before we summarize our main findings in Sec. |V[ 



II. MATERIALS AND METHODS 

A long-range attractive interaction was induced be- 
tween emulsion colloids by polymer depletion. The emul- 
sion consisted of poly(dimethyl siloxane) drops dispersed 
in a solvent mixture of 1,2-ethane diol (ED) and wa- 
ter (mass fraction of ED = 0.59). The solvent compo- 
sition was adjusted to closely match the refractive in- 
dex of the emulsion to minimize van der Waals attrac- 
tions between drops and to enable confocal imaging to be 
conducted deep within the sample. A particle radius of 
a = 316 ±11 nm and a size polydispersity of 0.17 ± 0.07 
was determined from dynamic light scattering measure- 
ments. The thickness of the polymeric stabilizing layer 
surrounding each emulsion drop was evaluated by cen- 
trifuging a suspension and equating the packing fraction 
of the sediment to the jamming density of a hard sphere 
system with the same polydispersity [20] . This procedure 
gave a layer thickness of « 7 ± 1 nm. The density mis- 
match between emulsion drops and the continuous phase 
is Ap = —130 ± 10 kg m -3 . To induce a depletion in- 
teraction, we added the non-adsorbing anionic polymer 
xanthan (Kelco, M w = 4.66 x 10 6 g mol -1 ). The polymer 
radius of gyration was determined as r g = 194±10 nm by 
light scattering and viscometry. The strength of the de- 
pletion attraction generated is a function of the polymer 
concentration and its range is controlled by the relative 
size r g /a = 0.62 ± 0.04 of the polymer and particle. The 
polymer concentration is quoted here in terms of the di- 
mensionless ratio c p /c*, where c* = 3M w /47rrg7VA is the 
overlap concentration (c* = 0.25 mg/ml) and N\ is Avo- 
gadro's constant. Full details of sample preparation are 
contained in Teece et al. [T7] . 

To monitor the collapse of the gels we used time-lapse 
video recording to record images of the emulsions as they 
cream. A low magnification image of the settling gel was 
projected onto a CCD camera (Allied Vision Technolo- 
gies F-080B). A regular sequence of images was captured 
every 20 seconds. The image series was corrected for op- 
tical distortion and non-uniformities in illumination be- 
fore being calibrated using an accurate grid of lines. The 
images near the center of the cells were analyzed and 
the interface separating the upper (dark) phase from the 
lower (bright) phase was identified automatically using 
an image analysis routine. The height h of the interface 
was extracted as a function of time with an accuracy of 
about ±0.3 mm. To aid visualization a low concentra- 
tion (« 0.001 mg/ml) of an adsorbing black dye, Sudan 
black, was added. The dye preferentially partitions into 



the index-matched PDMS drops so that the colloid-rich 
phase appears dark in transmitted light. The colloid- 
polymer mixtures were thoroughly mixed at the start 
of the experiments before being loaded into cylindrical 
glass vials with an internal diameter of d — 17 mm. The 
cell diameter was varied between 15-23 mm and both 
cylindrical glass and poly(styrene) cells were used, with 
no significant change in collapse behavior. To eliminate 
air bubbles which lead to irreproducible settling dynam- 
ics we used a gentle slow tumbling of the sample vial to 
thoroughly mix the samples before observation. Repeat 
experiments showed that following this protocol the col- 
lapse kinetics could be measured with a reproducibility 
of about 10-15%. 

Rheological measurements were performed at 23° C 
with a Bohlin HR Nano rheometer (Malvern Instru- 
ments). To study simultaneously the temporal evolution 
of the elastic properties and the height of the gel a novel 
rheometric vane experiment was developed which allowed 
visual observation of the gel as rheological measurements 
were performed. The vane was made from stainless steel 
and consisted of four blades (diameter 22.7 mm, height 
10 mm). The vane was carefully inserted into clear poly- 
carbonate sample vials (diameter 25 mm, height 65 mm), 
10 mm below the top surface of the gel and a thin layer 
of silicone oil added to minimize evaporation. The vane 
remained inside the dense upper phase for the duration of 
the rheological experiments, allowing the collapse process 
to be continuously monitored. Oscillation measurements 
were performed at 0.5 Hz at intervals of 200-250 s under 
controlled stress conditions, within the experimentally- 
determined linear viscoelastic region, while the height of 
the gel was monitored simultaneously by time-lapse video 
microscopy. The absence of wall slip was confirmed by 
watching the movement of small air bubbles deliberately 
introduced into samples. 

To directly probe changes in the microscopic topology 
of the gel during collapse we used fluorescent confocal 
microscopy. The continuous phase of the gel was la- 
beled with 0.02 mg/ml of the fluorescent dye rhodamine- 
B which combined with the high transparency of the 
emulsions provided by refractive-index matching allowed 
high resolution optical visualization deep within the gel. 
A light microscope (Zeiss, Axioskop S100) was mounted 
horizontally on its side, at right angles to gravity, and 
two-dimensional fluorescent images of regions 146 x 146 
fim were acquired at 543 nm. The gel was contained in 
a square cross-section glass vial, with an internal dimen- 
sion of 13 mm, mounted on a low profile translation stage 
so that the gel could be imaged at different vertical po- 
sitions, throughout the full 30 mm height of the sample 
[T7] . Since the emulsion drop radius is below the optical 
resolution limit of the confocal microscope we can not 
identify individual drops. Instead, we concentrated on 
the larger scale structure of the gel. The bicontinuous 
network was identified by thresholding the confocal im- 
ages to determine the location of the interface separating 
the (dark) emulsion phase from the (bright) continuous 
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FIG. 1. (Color online). State diagram showing gels stud- 
ied. The solid curve corresponds to the gas- liquid binodal 
calculated from the generalized free volume theory [22], for 
a polymer-colloid size ratio r g /a = 0.62. The open trian- 
gles identify the experimentally determined phase boundary. 
The theoretical prediction for the critical point is shown by 
the filled circle. The region below the coexistence line can 
be separated into two kinetic regimes: a region of complete 
demixing (I), and gelation (II). Quenched into region II, sus- 
pensions form a space-spanning network consisting of thick 
strands of particles. The vertical line indicates the constant 
colloid volume fraction gels studied (0 = 0.213). The colored 
symbols represent estimates of <j> for the strands of particles 
(blue squares) and the coexisting gas (red squares) after phase 
separation is complete locally. Error bars represent the age- 
dependent variation in from the same sample. 



phase. To correct for in-plane variation in the fluores- 
cence yield, each image was divided into 16 sub-images 
and a local threshold for each sub-image was determined 
using a cluster-based algorithm [21]. A careful analy- 
sis of the resulting binary images, backed up by direct 
observation, showed that this approach reliably located 
the shape and positions of the emulsion and aqueous do- 
mains. 



III. EXPERIMENTAL RESULTS 

A. Collapse dynamics 

To begin, the phase behavior of mixtures of emulsion 
and polymer was investigated as a function of both the 
emulsion volume fraction and polymer concentration 
Cp/c*. The state diagram plotted in Fig. [l] summarizes 
the results and shows the locations of a stable liquid 
phase, a narrow region of equilibrium gas-liquid demix- 
ing (I) and a broad zone of non-equilibrium gelation (II). 
The generalized free volume predictions (GFVT) [22] for 
a polymer-colloid size ratio r g /a = 0.62 in a good solvent 
are shown by the solid lines in Fig. [I] The agreement be- 
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FIG. 2. (Color online). Sudden collapse of a gel. (a) Time- 
lapse images of an emulsion-polymer mixture, with compo- 
sition (j) = 0.21 and c p /c* = 3.6, as a function of the time 
after shaking. Each image is labelled by the time elapsed t w , 
in units of 10 3 s. The initial height ho of the sample is 40 
mm. The characteristic delay time after which the net- 
work starts to collapse is indicated by the dashed arrow. The 
solid line denotes the position of the interface between the 
upper colloid-rich and lower polymer-rich phases, (b) The 
normalized height h/ho of the gel shown in part (a) as a func- 
tion of the elapsed time showing the three stages of settling 
characteristic of sudden collapse. 



tween the calculated gas-liquid binodal and experiments 
is good confirming that the experimental system is accu- 
rately represented by a simple mixture of hard spheres 
and non-adsorbing polymer chains. 

The sudden collapse of gels was investigated as a func- 
tion of both the strength of the attractions —U c /k^T and 
the initial height ho of the sample. In the absence of poly- 
mer, emulsions remained stable and showed no noticeable 
separation so the mechanical instability seen is a conse- 
quence of aggregate formation. The process of collapse is 
exemplified by the time-lapse CCD images reproduced in 
Fig. [2|a). Qualitatively we identify three distinct stages, 
characterized by the interface velocity v = d/i/dt w , where 
t w is the age of the gel. During an initial lag period (I) 
the network of attractive particles produces a mechani- 
cally stable solid, which is capable of supporting its own 
weight. This regime of solid-like stability persists how- 
ever only for a limited duration. On times longer than 
Td, the lag time, the network yields and a clear interface 
appears (identified by the solid line in Fig. |2^a)). The 
interface velocity v grows smoothly as the gel shrinks 
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and the collapse becomes progressively more rapid. This 
period (II) of rapid collapse terminates when phase sep- 
aration nears completion and the interface approaches 
the final equilibrium plateau. In the final consolidation 
stage (III), the settling velocity drops markedly as the 
collapsed gel continues to slowly compress like a solid 
under its own weight. 

The lag time Td is a strong function of the polymer 
concentration and hence the strength of depletion attrac- 
tions. Fig [3^a) shows the time-dependent height h(t w ) 
of gels prepared with different polymer concentrations 
but for the fixed initial height ho = 40 mm. Inspection 
reveals two striking features. First, as reported in pre- 
vious work [10] [TT], Td grows strongly with increasing 
polymer concentration. Indeed the concentration depen- 
dence of the lag time is well described by the exponen- 
tial relationship, Td ~ exp(c p /c*), as shown in Fig. [3|c). 
Second, the sedimentation profiles are remarkably simi- 
lar in shape when plotted in a linear-log representation. 
The height profile at a low polymer concentration may 
be simply mapped onto a high concentration sample by 
shifting the collapse profile to the right along the loga- 
rithmic time-axis. To explore this scaling behavior quan- 
titatively, we focus purely on the collapse regime and 
replot the change Ah = ho — h in the height of the gel 
as a function of the shifted time variable r = t w — Td, 
the time counted from the instant when the network first 
yields. Fig.|3jb) demonstrates, rather unexpectedly, that 
the initial change in the height of the gel depends linearly 
upon t 3 / 2 , over a wide range of polymer concentrations. 
In the early stages of gel collapse, where Ah < 0.5/io, 
the height of the gel follows approximately the algebraic 
expression 
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h(r) = h - Ari 
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where A is a polymer-dependent prefactor. The exis- 
tence of this simple relationship suggests a single com- 
mon mechanism is controlling the collapse of the gel at 
different polymer concentrations. We return to a more 



detailed discussion of this point in Sec. [IV] to speculate 
on a possible origin of the r 2 scaling found here. Finally, 
we note that an alternative scaling has been suggested 
by Kilfoil et al. [23 . We found however that their ap- 
proach failed when applied to the wide range of polymer 
concentrations and heights studied here. 

Previous studies [13 , 24 of gel collapse have suggested 
that the mechanism of collapse depends sensitively on the 
initial height ho of the gel. Gels formed in short sample 
cells display steady or 'creeping' sedimentation where the 
height falls continuously with age at a rate which decays 
exponentially with time while taller samples show sud- 
den collapse. To test whether this behavior is intrinsic 
to the long-range systems studied here we have varied 
the initial height ho and monitored the evolution of the 
height of the gel with time. In all the samples reported 
here, polymer concentration from c p /c* = 2.4 — 4.0 and 
heights ho = 22 — 63 mm, sudden gel collapse was always 
observed and we saw no transition to creeping sedimen- 
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FIG. 3. (Color online). Gel collapse at fixed height, (a) Time 
evolution of gel height for different polymer concentrations. 
Only data in the lag and collapse regimes is shown for clar- 
ity, (b) The change Ah in height as a function of t 3//2 , where 
r = t w — Td is the time elapsed from the start of collapse. 
Symbols and data are the same as in (a). The solid lines are 
straight line fits. The inset shows an expanded view of the 
short-time data, (c) Comparison between the experimentally- 
measured lag times Td (filled points) and the average lifetime 
Tesc of an individual particle bond (solid line). Lag times 
were measured in both glass (triangles) and poly(styrene) 
vials (squares). The dashed line, which reproduces the ex- 
perimental data reasonably well, equates to a fixed number 
of bond lifetimes (t& psj 240 r eS c)- 
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FIG. 4. (Color online). Height-independent lag times, (a) 
Temporal evolution of a gel with different initial heights but 
a fixed composition, <j) — 0.213, and c p /c* = 3.6. The arrows 
indicate the lag time r^. (b) The lag time as a function 
of the initial height ho of the gel, for a number of different 
polymer concentrations. The open symbols denote data ob- 
tained in glass-walled cells while the filled symbols indicate 
measurements in poly(styrene) cells. The nature of the cell 
wall has no noticeable effect on the delay time measured. 



tation. Fig.[4^a) shows a representative set of data where 
the time evolution of the interface height h(t w ) is plotted 
for a range of initial heights and a single polymer con- 
centration (cp/c* = 3.6). Inspection of the data reveals 
that rather surprisingly the time during which the gel 
is solid-like is largely independent of the initial height 
ho of the sample. Measurement of the variation of ra 
with height for a wide range of polymer concentrations, 
presented in Fig. Qb), confirms this observation. We 
have checked that this is not due to solid friction between 
the gel and cell wall [19] by repeating measurements at 
Cp/c* = 2.4 using cylindrical poly(styrene) cells to alter 
the degree of wall adhesion. The results, shown as the 
filled data points in Fig. [4jb), are in excellent agreement 
with the data obtained in the glass vials (open points), 
demonstrating that wall friction is unimportant. 

Since the gel is initially a solid, the top of the sam- 
ple vial is subject to a gravitational stress <r g , which is 
generated by the full buoyant stress of all of the sus- 



pension below, so a g = Apgcf)ho. Taking <p = 0.21 and 
Ap = —130 kg m -3 we estimate that a g varies from be- 
tween approximately 5 to 15 Pa for the heights used here. 
The values for the buoyant stress considerably exceeds 
the yield stress of the gel network, which we estimate 
from rheological measurements as a y ~ 0.1 Pa, so even 
while cr g is an order of magnitude larger than the stress 
required to break the gel's load-bearing structure the gel 
does not collapse macroscopically. The insensitivity of 
Td to ho means we can rule out the possibility that the 
initial period of latency of the gel is determined purely 
by the breaking of single uncorrelated bonds. Collapse 
clearly requires a substantially larger degree of restruc- 
turing of the network than is necessary simply for me- 
chanical yielding. 

To explore the effect of height on the kinetics of col- 
lapse we focus on the initial rate of collapse of the gel 
with different ho- Figure [5] shows that the power-law ex- 
pression (Eq. [TJ which captures well the height variation 
in gels with a fixed ho also holds for gels with a wide 
variety of different starting heights. The gel does not 
collapse with a fixed time- invariant velocity but rather 
the interface velocity v = dh/dr behaves at short times 
like t 1 / 2 , a behavior which hints at a surprisingly novel 
mechanism of collapse. To further investigate this mech- 
anism we have studied the dependence of the prefactor 
A = dA/i/dr 3 / 2 on the height and polymer content of the 
gel. Fig. pfc) shows that rather remarkably, for all the 
gels studied, the prefactor A does not change with the ini- 
tial height of the gel. Since the gravitational stress a g on 
the gel increases linearly with its height, this result sug- 
gests that <7 g is relatively unimportant in determining the 
process of collapse, at least under the conditions of our 
experiments. The central role of thermally-induced bond 
dissociation is seen in Fig. [5jd) where it is revealed that 
the prefactor A scales exponentially with c p , equivalent 
to an exponential dependence on the depth of the inter- 
action potential. In conclusion, the process of collapse 
appears to be thermally rather than stress-activated. 

To interpret these striking observations we model the 
initial deformation of a gel using the poroelastic formal- 
ism first introduced by Buscall and White [25]. In this 
approach the gel is treated as a biphasic fluid-saturated 
porous continuum with the pore pressure P as a state 
variable. A gel consists of two distinct phases: a solid 
phase of connected strands of emulsion particles, and a 
second liquid phase consisting of a fluid solution of a non- 
adsorbing polymer. In the early stages of collapse the gel 
is essentially uncompressed so there is no elastic stress 
due to deformation and the rate of collapse is limited es- 
sentially by the rate at which fluid is forced out of the 
gel [18]. Defining v as the macroscopic velocity of the 
fluid flow through the gel and w as the local displace- 
ment of the solid network along the gravitational z-axis 
then using Darcy's law, 
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where k is the permeability of the network, which since 
we are considering only the early stage of collapse we 
assume to be height independent, and 77 is the viscosity 
of the continuous phase. Continuity demands 

(1 - <t>)v = -<\>d t w (3) 

which, since <p < 1 — </>, implies that the fluid velocity \v\ 
must be small [18 in comparison to \dtw\. Consequently, 
if the displacement of the gel varies as r 3//2 then, from 
Eq. [2j the pressure gradient at the top of the gel, adja- 
cent to the interface, must be increasing as r 1 / 2 . This 
time dependence rules out a simple compression of the 
gel as a consequence of gravity because the pressure gra- 
dient would then be a constant, d z P = —Apgcj), and the 
gel would accordingly shrink linearly with time [18] . The 
r 1 / 2 dependence of P suggests instead a diffusive process 
may be responsible for the unusual collapse dynamics 
seen. We can however rule out a bulk diffusive process of 
conventional syneresis, akin to the shrinkage of a poly- 
mer gel undergoing a phase transition [26 . In this case, 
the contraction of the matrix as the phase separation 
ensues would generate an increase in P which leads to 
the expulsion of fluid and a shrinkage of the gel as t 1 / 2 
[27] rather than the t 3 ^ 2 dependence seen here. Further- 
more, this process would require fluid to be transported 
through the full length of the system so the rate of diffu- 
sion would depend on the height of the gel, which is also 
incompatible with our data. Clearly, a new mechanism 
is required to correctly explain the observed data. 

B. Microstructure 

Sudden collapse reflects a dramatic loss of mechanical 
integrity as the gel ages. To understand the nature of this 
mechanical failure, we first examine the evolution of the 
elastic properties of the gel while simultaneously record- 
ing the height of the gel. Figure [6] shows the rheology 
and interface height h(t w ) during the settling of a repre- 
sentative gel. At the earliest times recorded, the sample 
is solid-like with an elastic shear modulus G' which is 
larger than the viscous modulus G" . The measurements 
however reveal that the elasticity of the gel far from re- 
ducing with time, as one might naively expect, actually 
increases continuously up to the point t w = when the 
gel starts to collapse. Immediately collapse starts, G' also 
drops, but only by a relatively small amount (less than 
10%), before finally growing steadily as the emulsion vol- 
ume fraction in the upper phase increases with the onset 
of phase separation. We see no sign of large-scale hy- 
drodynamic mixing, recirculation, and the development 
of channels which have been seen in some other stud- 
ies of gravitational collapse [8] [12, 13 . It is intriguing 
that even though G' drops at the initiation of collapse 
the overall mechanical response remains predominantly 
elastic (G' > G"). We suspect that this reflects the in- 
homogeneous nature of the network at collapse. If the 
sample is heterogeneous the response measured will be 



an average over regions where the gel has broken apart, 
and is a fluid, with other portions of the network which 
remain elastic. The recorded response will then depend 
on the relative sizes of the mechanical vane and the in- 
homogeneous regions within the gel. To confirm that the 
sample becomes a fluid when the gel begins to collapse 
we have placed a small glass block about half the width 
of the cell in the base of the cell. When the gel begins 
to collapse the interface between the gel and the bottom 
of the cell rapidly flattens indicating that the base of the 
gel becomes a fluid of aggregates as the gel collapses. 

While the rheological measurements provide a mechan- 
ical insight into gel settling, they do not clarify the link 
between the macroscopic processes of collapse and the 
microscopic structural reorganization occurring during 
aging and sedimentation. Indeed, at first sight, it seems 
counter-intuitive that a gel which is becoming gradually 
stiffer with time should ever collapse at all. To probe 
the link between the microscopic and macroscopic length 
scales, we have examined the temporal evolution of the 
gel microstructure using confocal microscopy. The coars- 
ening is illustrated by the binary two-dimensional im- 
ages reproduced in Fig. [7] In the rectangular cell used 
for imaging experiments, the delay time was measured 
as Td ~ 1.3 x 10 5 s so both images refer to the latency 
period before collapse starts. Clearly, although the gel 
remains mechanically stable during this period there is 
a slow but continuous evolution in the microscopic na- 
ture of the particle network and the system is not ar- 
rested. A closer look at the data in Fig [7] reveals that 
the interface between the continuous and particle phases 
is quite rough, suggesting that surface tension is unim- 
portant and the dense portion of the gel is not a fluid. 
Direct observation show that particle diffusion is strongly 
suppressed and particles move infrequently between the 
strands of the network, indicating that the interaction 
network is probably glassy. The images in Fig [7] illustrate 
two further microstructural characteristics which will be 
important to our later discussion on the mechanism of 
network collapse. First, it is evident that the thickness 
of the backbone of the gel grows slowly but continuously 
with time, a fact which probably explains the increase 
of G' with t w seen in Fig. [6j Second, the thickness of 
the network of particle strands is not uniform. The gel 
contains a relatively large number of thin junction points 
where two or three arms (in 2D) are connected together. 
Simultaneous breakage of the particle bonds at these rel- 
atively weak junction points would lead to a rapid break 
up of the whole network. 

To quantify the change of the microstructure with 
time, we use chord methods developed to analyze sta- 
tistically random heterogeneous materials [28 . We su- 
perimpose on the images of Fig [7] a uniform grid of hori- 
zontal and vertical lines. The two-phase interfaces divide 
each grid line into chords that are either inside the dense 
part of the network or else lie within the solvent back- 
ground. We define a chord as the line segment between 
two consecutive intersections of the random line with the 
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FIG. 5. (Color online). Initial r 3/ ^ 2 -collapse dynamics, (a) The height h of a gel plotted as a function of the 3/2-power of 
the time elapsed after the gel yields, for different initial heights ho. Curves are labelled by the initial height. The samples 
have fixed polymer concentration c p /c* = 3.6, and colloid content <j) — 0.213. (b) Similar time-dependent settling observed in 
gels with Cp/c* = 2.4. (c) Invariance of the prefactor A — — lim T ^o dh/dr 3 ^ 2 with the height of the gel. Curves are labelled 
by Cp/c*. (d) Exponential dependence of A on the polymer concentration, highlighting the activated nature of the collapse 
process. 



bounding two-phase interface. Focusing only on those 
chords that lie within the dense strands of the network, 
we count the number of chords N(l) with lengths in the 
range I and I + dl. If TV is the total number of chords 
then the degree of linear 'connectedness' of the gel may 
be characterized in terms of a probability density func- 
tion, p(Z), where 



p(l)dl 



N ' 



(4) 



The quantity p(l)dl is the probability that a randomly- 
chosen chord has a length between I and I + dl. We find 
that the chord length distribution p(l) displays a char- 
acteristic shape, with p(l) first increasing with growing I 
before reaching a maximum at a finite I and finally de- 
caying exponentially for larger I. The diameter of the 
strands of particles within the gel may be characterized 
either from the value at which p(l) takes a maximum 
value or from the mean chord length l g 



I* 



poo 

/ lp(l)dl, 
Jo 



(5) 



which we use here because it displays a smaller statistical 
error. Finally, the chord functions also provide an effi- 
cient means to estimate the volume fraction of the high 
density colloidal regions in the gel. If we assume the gel 
is isotropic and the two-dimensional images are chosen 
randomly then the fraction \ °f the volume of the gel 
occupied by the dense regions is 



X : 
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(6) 



where L is the total length of the original lines. The 
volume fraction of colloids, g , in the dense regions of 
the gel is given by the ratio </>/x? where <\> is the initial 
colloid volume fraction. 

Next we characterize the slow evolution of the gel 
structure. Fig. [8] shows the age dependence of the mean 
chord length l g and the volume fraction g prior to col- 
lapse. A key observation is that during the latency period 
when the height of the gel is unchanged structural reor- 
ganization is never fully arrested but continues, albeit 
rather slowly. So, for instance, the growth of the mean 
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FIG. 6. (Color online). Simultaneous measurements of lin- 
ear viscoelasticity and height h of gel as a function of time 
elapsed since preparation t w . Gel had an initial composition 
of Cp/c* = 3.6, cj) — 0.213 and a measured delay time of 
= 1.2 x 10 5 s (indicated by dashed line). The elastic G' 
and loss G" moduli were measured by applying a oscillatory 
stress of magnitude 0.0025 Pa at a frequency of 0.5 Hz and 
measuring the strain response. The gel stiffens continuously 
with age until the elastic modulus drops at t w = as the gel 
begins to collapse. 




FIG. 7. Continuous coarsening of gel with age. Two- 
dimensional binary representation of a gel with composition 
c p / C ; = 2.4, (j) = 0.213 after (a) 3.6 xlO 3 s, and (b) 1.12 x 10 5 s 
from preparation. Gravity points vertically downwards and 
the scale bar corresponds to 20 /mi. 
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FIG. 8. (Color online). Evolution of colloidal network with 
time, (a) Average chord length l g in the colloid phase, in units 
of the particle diameter. The rate of increase of / g is always 
slower than the diffusive growth, l g ~ t\l 3 (shown by dashed 
line) characteristic of classical spinodal decomposition [29] , 
and slows down considerably with increasing polymer concen- 
tration Cp/c*. (b) Average colloid volume fraction g within 
dense regions of the gel, as a function of elapsed time. The 
density of the particle strands within the network increases 
progressively with age before rising rapidly at the onset of gel 
collapse. 



C. Origin of lag time 



chord length is well described by a power law, Z g ~ t^, 
with an exponent a. The growth law is always much 

1 /3 

slower than the tj dependence expected for the diffu- 
sive regime of classical liquid-gas phase separation. The 
growth exponent a depend rather strongly on the poly- 
mer concentration with a decreasing markedly as c p /c* 
is increased. Similar slow growth has been identified pre- 
viously in simulations of deep quenched Lennard- Jones 
systems [29] and has been interpreted as indicating that 
the dense domains of the gel are actually glassy. 



Having characterized the macroscopic process of col- 
lapse, we now discuss the mechanism for the initial failure 
of the particle network. A gel is a metastable phase with 
a high free energy density whose consolidation is driven 
ultimately by the thermodynamic driving force for phase 
separation. However once a stable percolating network of 
strands of particles is formed the dynamics of phase sep- 
aration is slowed down considerably because, as evident 
from Fig. [7| the strands of the network are many parti- 
cles wide so large scale rearrangements of the gel require 
the simultaneous dissociation of all of the particle bonds 
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within the cross-section of a strand, which will be very 
rare. The network accordingly lowers its free energy via 
a series of small-scale structural reorganizations which 
proceed through the rupture of essentially single parti- 
cle bonds. The breakup of an energetic bond between 
particles, diffusion to dense region of the network, and 
a reformation of the broken bond allows a net increase 
in the number of nearest neigbouring particles with a 
concomitant lowering of the free energy of the system. 
For the network to coarsen, the system must overcome 
the energetic barrier associated with single bond rupture. 
This could be achieved, in principle, either thermally or 
as a result of an applied stress. The observation that 
the delay time is unaffected by the initial height of the 
gel strongly suggests that the delay time and hence the 
rupture of individual bonds is controlled primarily by 
thermal fluctuations rather than being stress-driven. 

To calculate the average lifetime r esc of an individ- 
ual particle bond due to thermal fluctuations, we assume 
that a single bond ruptures on a scale comparable to 
the time it takes a Brownian particle to escape from an 
attractive ramp potential with the same range S/a and 
depth — U c /k^T as the interparticle depletion potential. 
The mean escape time in the overdamped limit is given 
by the Kramers expression [30] 

_ S 2 exp(-U c /k B T) - (1 - U c /k B T) 

Tesc A {u c /k B Ty [n 

where D t is a translational diffusion constant. We esti- 
mate D t from the short-time self diffusion constant in 
a hard sphere suspension at the same 0, which since 
the dense regions of the gel have a volume fraction 
~ 0.55 is about 20% of the dilute free particle limit, 
Do = ksT / (6iri%a) . The limiting low shear viscos- 
ity ?7l was determined by extrapolating measurements of 
the steady-shear rheology of the polymer solution to a 
vanishing shear rate and fitting to the Martin equation, 

— = 1 + [r]}c p exp (fc H Mcp) (8) 
Vo 

which has been found to correlate well viscosity in dilute 
and semi-dilute polymer solutions (c p /c* < 10). Here [rj] 
is the intrinsic viscosity, r]o is the viscosity of the mixed 
solvent, c p the polymer mass concentration, and &h is 
a constant (equivalent to the Huggins constant at low 
c p ). Fitting data in the range c p = 0.6 — 1.2 mg/ml 
to Eq. [8] gave [rj] = 2.32 ml/mg and fc H = 1.2. The 
width of the depletion zone S and the potential at con- 
tact — U c /k^T were estimated using the generalized free 
volume theory for mixtures of hard sphere colloids and 
excluded- volume polymer chains, as detailed in Ref. [T7] . 
Figure [3]^c) shows a comparison between the measured 
delay time and the average lifetime r esc of a single 
particle bond, estimated from Eq. [7J The ratio of the 
two timescales is very nearly constant, for a wide range 
of polymer concentrations, with the delay time approxi- 
mately 240 times the estimated single particle Kramers 



escape time. This strong correlation highlights the piv- 
otal role of spontaneous thermal fluctuations in determin- 
ing the latency period of the gel. The fact that the delay 
time is many times longer than the rupture of a single 
bond probably reflects the cooperative nature of gel fail- 
ure. The strands of the network are several particles wide 
so failure requires the simultaneous dissociation of all of 
the bonds in the cross-section of a particle chain [31]. 
The alternative picture proposed by Buscall et al. [32] , 
that the ratio Td/r esc is determined by the mean coordi- 
nation number of particles within the gel, could only be 
consistent with our observations if the mean particle co- 
ordination number varied with the depth of the attractive 
potential. To distinguish completely between these two 
possibilities requires a more detailed microscopic model 
of gel failure than is currently available. 

To explore the microscopic mechanism by which ther- 
mal fluctuations lead to delayed failure we used real-space 
confocal imaging to follow the time evolution of the gel. 
Since the load-bearing nature of the network is clearly 
important we concentrate on changes in the connectiv- 
ity of the strands of particles which constitute the gel. 
Figure |9ja) shows two-dimensional confocal slices taken 
from the same physical region within an aging gel be- 
fore collapse occurs. Because of the finite bond energy, 
the network structure slowly but continuously evolves, 
with fluctuations in both the number and type of junc- 
tion points. By comparing 2D images of the fine-stranded 
structure of the network at hourly intervals we identi- 
fied discrete strand association and dissociation events 
occurring over this period. Examples where the strand 
network is ruptured are indicated by the dashed circles in 
Fig. [9^ a) while the solid circles identify new cross-links 
formed by the reassociation of strands. Counting the 
number of reassociation AT link and rupture iVb re ak events 
recorded per hour, as a function of the age of the gel, 
results in the data shown in Fig. [9^b). There is gradual 
reduction over time in the number of reassociation events 
Mink? as the network is formed in an open high-energy 
state and then relaxes slowly into a lower more com- 
pact structure. Strikingly however, we see that the rate 
of bond rupture does not show the same slowing-down. 
^break is essentially independent of age, presumably be- 
cause rupture is an activated process which is dominated 
by the single particle bond energy barrier. The conse- 
quence of the different time dependence seen for asso- 
ciation and rupture is that the proportion of breakage 
events / = iV break /(AT break + JViink) (shown in the inset of 
Fig. [9jb)) grows with the age of the gel. The increasing 
proportion of strand ruptures ultimately leads to failure 
of the stress-bearing backbone of the gel and the initia- 
tion of collapse. 



D. Appearance of structural heterogeneities 

Work in the last decade [33j [34] has shown that 
soft glassy materials frequently display structural het- 
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FIG. 9. (Color online), (a) 2D confocal images of coarsening gel network formed at <j> — 0.213, and c p /c* = 2.4. The particles 
are shown in white. The solid circles indicate network junctions which have formed in the preceding 3600 s. The dashed circles 
indicate positions where the network has, in the same interval, broken. The scale bar is 30 fim long. Network collapse occurs 
at Td ~ 120 x 10 3 s. (b) The number of reassociation and rupture iVb rea k events per hour as a function of the age t w of 

the gel. The dashed lines are guides to the eye. The inset shows the t w -dependence of the fraction / of rupture events. 



erogeneities. In materials where the elastic behavior of 
a material dominates over its viscous response any de- 
formation due to a local rearrangement can propagate 
macroscopic distances so the size of regions which un- 
dergo correlated rearrangements can be sizeable. If this 
holds true in our system, then it should be feasible to see 
signs of the long-range distortion field generated by local 
rearrangement events by microscopy. 

To test these ideas we have used confocal microscopy 
to monitor the time evolution of the network structure 
as a function of the vertical z-position within a gel. A 
series of 2D-confocal images were collected at regularly- 
spaced 1 mm heights from a colloid-polymer gel with a 
total height of ho = 15 mm. Images were acquired for 
~ 7 hours after the cessation of mixing, until the point at 
which gel collapse occurred. The characteristic domain 
size of the network R c (h,t w ) at a height h and time t w 



was calculated from the static structure factor S(q, t w ) 



S(q,t w ) 



I 



2nqAq Jq<\q>\< q+ Aq 



dq' (I(ci',t w )I{-Ci',t w ) 



(9) 

where /(q, t w ) is the two-dimensional Fourier-transform 
of the image intensity J(r,t w ) at time t w , Aq = 2tt/W, 
and W is the image width. The domain size is R c = 
7r/ (q) where (q) = J dq qS(q)/ J dqS(q). Measure- 
ments of R c for different sample ages t w are plotted in 
Fig." 
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and confirm that the aging of the gel network 
shows considerable spatial diversity: the domain size is 
large in some regions of space and small in others. Im- 
mediately after mixing, we observe the formation of a 
very uniform network with an average domain size of 
(R c ) = 17.5 fim and a spatial variation of just 1.4% 
(standard deviation/ (R c ))- But after t w = 2 h, while 
the mean size has grown only slightly to (R c ) = 21.3 jam 
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FIG. 10. (Color online). The domain size R c in a gel of 
height ho — 15 mm (c p /c* = 2.4) as a function of the distance 
measured from the base of the cell. Note aging is hindered at 
the boundaries of the gel. The effect continues over a distance 
A which is of order a few millimeters. 



the spatial variation in R c has increased to 3%. After 
5 h, the variation in the domain size has increased still 
further to 9% ((Rc) = 26.7 fim). Clearly the data reveals 
that aging of the network is heterogeneous with spatial 
variation increasing with sample age. 



IV. DISCUSSION 

The most striking feature of our results is the appear- 
ance during collapse of the | -power law dependence of 
the height h(r) on the elapsed time r. A natural ques- 
tion is the physical origin of this unusual behavior. We 
propose that the dominant mechanism for collapse is the 
appearance of random micro-collapsed regions through- 
out the gel. On a microscopic level, the basic idea is 
that the particles comprising the gel attract each other 
relatively strongly so over time the gel spontaneously re- 
structures locally to create small regions of more dense 
packing. Since the collapsing particles are attached quite 
strongly to strands of the network as they rearrange they 
induce a local pressure field. This induced field, as a con- 
sequence of the poroelastic character of the gel, expands 
relatively slowly into the bulk of the gel. It is this long- 
range pressure field which we hypothesize generates the 
characteristic collapse dynamics evident in our experi- 
ments. Similar arguments have been invoked to account 
for the anomalous microscopic motion evident in dynam- 
ical light scattering of colloidal gels [33] [34] but not, as far 
as we are aware, for the macroscopic settling dynamics of 
gels. While our discussion has some features in common 
with the purely elastic models used previously [33 , 34 , we 
focus here on deformations at large (macroscopic) length 
scales where poroelastic fluid flows are important. 

For the moment, we idealize the gel as a one- 
dimensional chain of particles. Then if two particles 



leaves their equilibrium positions to stick together the 
left-hand particle imposes a force +/o on the left part 
of the chain at z' \ while the right-hand particle imposes 
an equal and opposite force — fo on the right-hand side 
of the chain located at z' + 5z' . The net effect of the 
rearrangement is therefore the creation of a local dipo- 
lar force f(z',t') at the random position z = z' . The 
intensity of this dipolar force, the dipole moment /i, is 
the product of the force /o and the displacement vector 
5z' in the limit as 5z' —> 0. At the dipolar stress center 
the fluid pressure P rises rapidly to a high value while 
further away P is almost unchanged. The resulting pres- 
sure gradient drives a flow through the porous medium 
and as fluid exits from around the applied force the gel 
deforms locally and more of the load is borne by the net- 
work. Eventually, the pressure at all points reaches the 
same value and the gel relaxes so that the applied load 
is everywhere balanced by the elastic stresses in the net- 
work. The time scale for this equilibration is determined 
by a diffusion constant D g1 with the deformation in the 
gel occurring over a length scale « y/D g t in a time t. 
The gel diffusion constant D g is [35] [36] 



Kk 



^(1-0)" 



(10) 



where k is the permeability of the network, rj^ is the vis- 
cosity of the continuous phase, and K is the bulk modu- 
lus of the network. Using measurements of the low shear 
viscosity (77L — 0.1 Pas), the elastic modulus of the gels 
(K ~ 10 Pa), and literature values for the permeability 
of similar density gels [18 (k ~ 10a 2 ~ 10 -12 m 2 ) we es- 
timate a diffusion constant in our system of D g ~ 10 -10 
m 2 s _1 . Both the pore pressure and, in general, the per- 
meability will change with time as the gel contracts lo- 
cally. However in the initial stages of collapse the gel 
is uncompressed, the permeability is constant, and the 
equations of linear poroelasticity apply [37]. We shall 
ignore all non-linear effects. The bulk shrinkage of an 
unconstrained gel is linear so it is natural to assume that 
the dipole intensity will also be a linear function of time, 
fi(i) fi t. 

For simplicity, we first consider the isotropic deforma- 
tion produced in the gel by the supposition of three con- 
tinuous orthogonal stress dipole centers (a single center 
of compression [38]). Rudnicki [39] has shown that the 
pore pressure at a distance r from a single continuous 
center of compression in a fluid-saturated porous solid is 
of the form 



P s (r,t) 



Q 



-erfc(£/2), 



(11) 



where £ = r/( J D g t) 1 / 2 , Q is proportional to /io, and erfc 
is the complementary error function [40 . As t — » 0, 
ercf(r )2^J D g t) —> 0, and the pressure is zero. At finite 
times, the fluid has an opportunity to diffuse away from 
the origin, and the spherically symmetric pressure wave 



expands. In Fig. 11 we plot a time series of the spread- 
ing pressure field as it diffuses away from the origin. The 
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length scale where the pressure is finite is controlled by 
fluid diffusion within the gel and is thus time-dependent. 
From the figure it is clear that the spatial extent of the 
pressure field is of order 6( J D g t) 1 / 2 . Finally, as t — >• oo, 
evcf(r/2^D g t) — > 1, and the pressure field approaches 
the pure 1/r-elastostatic solution, expected for a contin- 
uous dipole source [37] . 

We assume that at the start of collapse the gel contains 
centers of compression which are randomly distributed 
throughout the height of the gel, with a certain number 
p per unit length. To calculate the velocity of the gel in- 
terface dfW we note that, from Darcy's law (Eq.|2|, dtw is 
fixed by the total pressure gradient at the bottom of the 
gel. If we assume that each micro-collapse contributes in- 
dependently then the pressure gradient is simply the sum 
of the gradients generated by individual events occurring 
at different heights Zj within the gel. For simplicity we 
ignore the true vectorial nature of the problem and use a 
scalar approximation in which we assume that the pres- 
sure field due to each independent micro-collapse adds 
together coherently to produce the macroscopic pressure 
field. Since the vertical pressure gradient must control 
the process of collapse we focus solely on the z-direction. 
A single center of compression located a distance Zj from 



the gel interface creates a pressure gradient, Oj 
which from Eq. [TT] is 

-6j(zj,t) - — — 



d r P, 



^z)D, 



-erfc(- 



=) 



Qe 



4tt 3 / 2 ^£»; 



3/2 



(12) 



The growth data in Fig. [10] suggests that coarsening of 
the network is suppressed near a surface, i.e. micro- 
collapse events appear preferentially at distances z > A 
away from a surface. This seems plausible since the en- 
ergy barrier for a rearrangement near an interface will 
probably be higher than for the same event in the bulk 
of the gel, as the strain induced by the creation of the 
dipole is larger. Inspection of the data in Fig. [l0| suggests 
that A is of order a few millimeters. Hence, we assume 
that micro-collapses are uniformly distributed over the 
interval from z = X to z = ho — A. The total pressure 
gradient at the base of the gel is therefore, 



d z p 



base 



h —X 



0(z,t)dz. 



(13) 



Since experimentally we have observed that the collapse 
process does not vary with the total height ho of the 
gel the upper limit of the integral can be extended to 
z = oo without significant error. After inserting Eq. [l2j 
the resulting integral can be performed exactly with the 
result: 

pQ 



d z P{t) 



1/2 



Vie 



pQ\n\ 



(14) 



In the regime where t ^> A 2 /4D g , which from our esti- 
mates for D g and A equates to t ^> 10 3 s and holds for 
all but the shortest times studied, the expression for the 
pressure gradient simplifies to 



d z P(t) 



PQ 



t>AV4D g 4X2^3/2 D 



1/2 



y/i. (15) 



On timescales t « 10 3 s the gel has typically not col- 
lapsed to any significant degree (see for instance the 
data in Fig. [3| so the permeability of the gel is not sub- 
stantially changed from its initial value and D g is time- 
independent. 

The confocal data, presented in Sec. |III C] reveals that 
micro-collapses first appear in gels with an age t w of « 
so the time t available for the propagation of the pressure 
wave in Eq. [15] is r = t w — ra . By combining Eq. [2| [3J 



and [151 we obtain 



Aft(r) 



pkQ 



47/(1 - (j))W/ 2 D l 



1/2 



r 3/2 



(16) 



where we have assumed that the flow of fluid through 
the network determines the initial rate of collapse, and 
t ^> A 2 /4D g . This expression is in very good agreement 
with the experimental results where a similar exponent of 
3/2 has been found, thus providing convincing evidence 
for our simple phenomenological model. In addition, our 
arguments predict that the coefficient of r 3 / 2 , which we 
identify with the scale factor A in Eq. [I] should be a 
system constant, independent of the initial height of the 
gel. This agrees with the height-independence seen in the 
experimental data presented in Fig.[5jc). We expect that 
the formation of micro-collapses is thermally activated so 
their number density p will be of the form 



po exp 



AE\ 



(17) 



where AE is an energy barrier for rearrangement. Since 
AE/k^T will scale with the strength of the interparti- 
cle potential U c /k^T one expects that the scale factor A 
will depend exponentially on the interparticle potential, 
or equivalently the polymer concentration, in agreement 
with the experimental data plotted in Fig.[5|d). Finally, 
we note that our model displays no explicit dependence 
on gravity since we hypothesize that gel collapse is a 
consequence of irreversible aging of the particle network. 
Gravity simply dictates the direction of gel collapse and 
we believe that, in the current case, the gravitational 
stress on the network is not sufficiently large to signifi- 
cantly enhance the thermal relaxation of particle bonds. 
However, this conjecture still awaits a direct experimen- 
tal proof. A systematic investigation of the collapse of 
transient gels as a function of the gravitational stress, 
which could be achieved by for example changing the 
density mismatch Ap or by using micro-gravity condi- 
tions, would confirm this prediction. Unfortunately, no 
such data is currently available although we plan in the 
near future to start such measurements. 
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FIG. 11. (Color online). An illustration of the pore pressure field generated by a continuous compression center placed at 
the origin of a gel. Times and distances are scaled so the diffusion constant D g is unity. The ring of radius r — 6( J D g t) 1 ^ 2 
parameterizes the spatial extent of the induced pressure field. 



V. SUMMARY 



We have studied the gravitational collapse of a colloidal 
gel by a combination of confocal microscopy, time-lapse 
video imaging, and rheology focusing particularly on the 
effect of the initial height ho of the gel and the strength 
of attractions U c /k^T. The gels are made of emulsion 
drops suspended in a refractive index-matched mixture 
of ethylene glycol and water, with a high molecular poly- 
mer added to induce a weak long-range attraction. For 
all systems, the height h(t w ) of the gel shows a charac- 
teristic two-step decay as a function of age t w : for t w less 
than the lag time ra the system resists gravity and there 
is no significant deformation, but for t w > the gel 
abruptly yields and collapses. The change in the height 
Aft = ho—h(t w ) of the gel during collapse has a number of 
distinctive features. First, we find that the initial degree 
of settling is well described by the expression, Aft ~ r 3 / 2 , 
with r the time counted from the moment when collapse 
first starts. Second, both the process of collapse and the 
lag time are independent of the initial height of the 
gel. Microscopically, the gel consists of a network of in- 
terconnected strands of particles. Confocal microscopy 
reveals that continuous restructuring of this network oc- 



curs which, with increasing age, leads to the breaking of 
bonds between particle strands and a progressive weak- 
ening of the network. The subsequent reduction in the 
large scale connectivity of the network eventually triggers 
a macroscopic collapse. Measurement of the microscopic 
structure of the gel during settling show that the age- 
dependent changes in the network are not isotropic but 
are concentrated around large inhomogeneities within the 
sample. We hypothesize that the collapse of the gel is de- 
termined by the rate at which fluid can be expelled from 
the gel. A simple phenomenological model of fluid flow 
driven by the formation of random compression centers 
within the gel correctly accounts for the behavior exper- 
imentally observed. 
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